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Abstract. We have conducted a survey of X-ray sources in XMM-Newton observations of M31, 
examining their power density spectra (PDS) and spectral energy distributions (SEDs). Our 
automated source detection yielded 535 good X-ray sources; to date, we have studied 225 of 
them. In particular, we examined the PDS because low mass X-ray binaries (LMXBs) exhibit 
two distinctive types of PDS. At low accretion rates, the PDS is characterised by a broken power 
law, with the spectral index changing from ~0 to ~1 at some frequency in the range ~0.01-1 
Hz; we refer to such PDS as Type A. At higher accretion rates, the PDS is described by a simple 
power law; we call these PDS Type B. Of the 225 sources studied to date, 75 exhibit Type A 
variability, and are almost certainly LMXBs, while 6 show Type B but not Type A, and are likely 
LMXBs. Of these 81 candidate LMXBs, 71 are newly identified in this survey; furthermore, they 
are mostly found near the centre of M31. Furthermore, most of the X-ray population in the disc 
are associated with the spiral arms, making them likely high mass X-ray binaries (HMXBs). 
In general these HMXBs do not exhibit Type A variability, while many central X-ray sources 
(LMXBs) in the same luminosity range do. Hence the PDS may distinguish between LMXBs 
and HMXBs in this luminosity range. 
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1. Introduction 

We present results from a survey of X-ray point sources in XMM-Newton observations 
of t he Andromeda Galax y (M31). M31 is the nearest spiral galaxy to our own, at 760 
kpc Ijvan den Borah 2000), and its X-ray population is dominated by X-ray binaries. 

Our work differs from previous surveys of M31 in two important ways. Firstly, we made 
power density spectra (PDS) from combined EPIC, background subtracted lightcurves 
of each source. This was made possible for the first time by the unprecedented sensitivity 
of XMM-Newton, allowing us to detect variability in M31 X-ray sources on ^100 s 
timescales. Secondly, we obtained the 0.3-10 keV luminosity of each source from best fit 
models to their spectral energy distribution (SED); previous works have used an assumed 
model to convert from count rate to luminosity in all but the brightest few sources. 

We have examined four XMM-Newton observations of the central region of M31, along 
with one obse rvation each of the North 1, North 2, North 3, South 1 and South 2 
fields (see e.g. IPietsch et al.ll200rij) . We performed automated source detection on these 
regions, and found 535 so urces that were observed by all three EPIC detectors (see 
IShaw Greening et al.l [2005. for details). However, we have not yet studied these sources 
in a uniform manner. 
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Table 1. Classifications of the X-ray population of M31 by power density spectrum (PDS) and 
spectral energy distribution (SED). The SED is classified by the photon index of the power law 
component, F, or by the temperature of the thermal component, kTss. 



PDS 




SED 


Classification 


Type 


A 


r -1. 5-2.1 


NS or BH LMXB in low/hard state 


Type 


A 


r -2.4-3.0 


BH LMXB in very high/ 








intermediate (steep power law) state 


Type 


B 


BB dominated 


BH LMXB in high/soft state 






kT BB -0.7-2 keV 




Type 


B 


Not BB dominated 


NS LMXB in high state 



2. Combining PDS and SED information 

X-ray binaries exhibit many phenomena on time-scales spanning milliseconds to years. 
Low mass X-ray binaries (LMXBs) are particularly interesting because they exhibit very 
distinctive power density spectra that depend more on the accretion rate than on the 
nature of the primary. Both neutron star and black hole LMXBs exhibit PDS that are 
characterised by a broken power law at low accretion rates; the spectral index ch anges 
from —0 to —1 at a break fr equency in the range —0.01-1 Hz l(van der KlisHl993l . We 
call this Type A variability l)Barnard et al J 120041 . At higher accretion rates, the PDS 
is described by a simp le power law oyer th e some frequency range; we refer to this 
as Type B variability l|Barnard et al ] |2004l) . Van der Klis (1994) suggested that the 
transition from Type A to Type B variability occurs at some constant accretion rate, 
and suggested a transition at — 1% Eddington. We have empirical evidence that suggests 
that the transitio n from Type A to Type B v ariability occurs at —10% Eddington in the 
0.3-10 keV band l(Barnard et alJl2004 fopO.^ . Since the Eddington limit is proportional 
to the mass of the primary, we expect black hole LMXBs to exhibit Type A variability at 
higher luminosities than neutron star LMXBs. We assume a maximum neutron star mass 
of 3.1 Mq and classify as a possible black hole LMXB any X-ray source that exhibits 
Type A variability at a 0.3-10 keV luminosity >4xl0 37 erg s~ x . 

We also obtained EPIC-pn SEDs for each source in the 0.3-10 keV band, and an equiv- 
alent background SED. Response matrices and ancillary response files were generated for 
each source. The resulting SEDs were grouped for a minimum of 50 count bin -1 from 
sources with >500 counts, and for a minimum of 20 count bin -1 for SEDs with <500 
counts. We fitted several spectral models to each SED and obtained the 0.3-10 keV lu- 
minosity from the best fit. If the source+background spectrum contained <150 counts, 
and/or the source spectrum contained <20 counts, then the luminosity was estimated 
from the 0.3-10 keV pn intensity; we assumed an absorbed power law model with line of 
sight absorption equivalent to 10 21 H atom cm -2 and a photon index of 1.7. 

By combining PDS and SED information we were able to classify the X-ray population 
of M31 to an un precedented degree; Table [H summarises our c l assific ation scheme for 
LMXBs (see e.g. iMcClintock fc RemillarS l2004 iBarnard et al"ll2005h . We do not yet 
ha ve such a scheme for HMXBs . For luminosity functions of the different regions of M31, 
see lShaw Greening et al 1 il2005h . 



3. Results 

The initial survey c overed the four obse rvations of the 63 brightest sources in the 
central region of M31 IjBarnard et alJl2005jl : we have yet to study the —130 remaining 
sources. Similarly, we have only surveyed 36 out of 87 sources in South 1 region; this 



^100 new X-ray binaries in M31 



121 



TYPE A TYPE B 




Frequency (Hz) Frequency (Hz) 

Figure 1. Power density spectra from the 2002 January XMM-Newton observation of two 
sources in the central region of M31. The PDS are averaged over 96 intervals of 128 bins, lasting 
666 s, and grouped geometrically. 

survey also predates the automated source detection. However, we have studied all 83 
sources in North 1 and all 43 sources in North 2. We have not yet studied any of the 
88 sources in North 3 nor the 41 sources in South 2. All the remaining sources will be 
surveyed over the next few months. 

We have examined the fo ur observations of the 63 X-ray sources in the core in most 
detail ()Barnard et al.ll2005|) : we were able to classify 107 out of 252 PDS as either Type 
A, Type B or flat; flat PDS had no significant power above 2, which is the expected power 
for Poisson noise when using Leahy normalisation. Figure ^ shows examples of Type A 
and Type B PDS exhibited by M31 X-ray sources during the -60 ks XMM-Newton 
observation on 2002, January 6. Furthermore, we found that t he Type B and flat PDS 
were both consistent with high accretion rate LMXBs in M31 l|Barnard et alJl2005|) . 

In the core, 40 out of 63 sources exhibited Type A variability in at least one observation; 
these are almost certainly X-ray binaries. A further 6 sources exhibited Type B but not 
Type A variability, making them likely LMXBs. In the North 1 and North 2 fields, 25 out 
of 83 and 1 out of 43 sources respectively exhibit Type A variability; in the South 1 field, 
9 out of 36 sources exhibited Type A variability. Of the 40 sources that exhibited Type 
A variability in the core, 1 3 did so at luminositi es >4xl0 37 erg s _1 , and were classed as 
bl ack hole cand i dates (see iBarnard et alJl2005(l . 

IPietsch et"alll|2005(l surveyed these same XMM-Newton observations, finding 856 point 
sources. They catalogued 7 LMXBs and 9 likely LMXBs, as well as 27 certain and 10 
likely associations of X-ray sources with globular clusters. These 53 sources are all likely 
LMXBs. We have surveyed 28 of these sources to date, and find that 10 of them exhibit 
Type A variability. Hence, 71 of the 81 X-ray binary candidates were newly identified 
by this work. . In Fig. HI we present a mosaic of GALEX images of M31 in the far 
ultraviolet (see e.g iThilker et alJl2005l and |http://www.galex.caltec.edu| ), with the X- 
ray sources that we have already studied superposed. Black dots represent X-ray sources 
that exhibited Type A variability in at least one observation, while white dots represent 
sources that have not. We note that the central region is dominated by black dots; 63% 
of the surveyed sources exhibited Type A variability in the core. Meanwhile, 25-30% of 
sources in the North 1 and South 1 regions, and only 2% of sources in North 2, exhibited 
Type A variability (see Fig. |2J • 

The majority of sources in the North 1, North 2 and South 1 regions have 0.3-10 keV 
luminosities <10 37 erg s _1 . Almost all the core X-ray sources that we have surveyed 
in this regime with classified PDS show Type A variability; hence we would expect all 
LMXBs in the M31 disc to exhibit Type A variability also. However, we only detect Type 
A variability in ~20% of the disc sources. Furthermore, we see that the disc population 
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Figure 2. X-ray sources that we have already surveyed, superposed on a mosaic of GALEX 
FUV observations of M31. Black sources exhibit Type A variability, and are hence likely LMXBs. 
White sources do not exhibit Type A variability. Out of the 225 sources surveyed to date, 75 
exhibit Type A variability. 



of X-ray sources largely follows the spiral arms, suggesting that they are HMXBs. We 
therefore suggest that it may be possible to distinguish between LMXBs and HMXBs, 
based on their PDS at 0.3-10 keV luminosities <10 37 erg s _1 ; the HMXBs are not 
expected to exhibit Type A variability as they are not disc fed. We expect most of the 
~130 faint sources in the Core to exhibit Type A variability. 
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